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A tight binding parametrization of local spin density functional band theory is combined with 
a dynamical mean field treatment of correlations to obtain a theory of the magnetic transition 
temperature and optical conductivity and T — spinwave stiffness of a minimal model for the 
pseudocubic metallic CMR manganites such a La\-xSr x MnOz- The results indicate that previous 
estimates of T c obtained by one of us (Phys. Rev. B61 10738-49 (2000)) are in error, that in 
fact the materials are characterized by Hund's coupling J « 1.5eV, and that magnetic-order driven 
changes in the kinetic energy may not be the cause of the observed 'colossal' magnetoresistive and 
multiphase behavior in the manganites, raising questions about our present understanding of these 
materials. 

PACS numbers: 



I. INTRODUCTION 

Pseudocubic manganites such as La\- x Sr x MnOz are 
of continuing interest to physicists and materials scien- 
tists for many reasons, including the possibility that the 
materials are 'half metals', in other words, have ferromag- 
netic ground states with fully spin polarized conduction 
bands. The possibility of 'half metallicity' was under- 
stood in the 1950s by Wollan, Kohler, Zener, Goode- 
nough, and Anderson 1-4 , but remains controversial to- 
day, with some calculations 5 and some 6 (but not all 7 ) 
experiments suggesting the presence of minority carriers. 

In the pseudocubic manganites such as 
Lai— x Sr x Mn03 the possible half metallicity arises 
from correlation and crystal-field effects involving the 
Mn d-levels, which in a naive formal-valence picture 
contain 4 — x d-electrons (see e.g Ref 2 for further 
discussion). In brief, in these materials the nearly cubic 
crystal field splits the d-multiplet into a lower-lying 
tig symmetry triplet and a higher lying e g -symmetry 
doublet. Correlation effects cause the tig levels to 
contain three electrons with parallel spin, producing 
an electrically inert S — 3/2 'core spin', while the 
remaining 1 — x electrons go into e s -symmetry orbitals 
and can move throughout the crystal subject to a Hund's 
coupling J which favors configurations in which an e g 
electron on a given site has its spin parallel to the core 
spin on that site 8 Evidently, if the Hund's coupling is 
strong enough, then a state in which all of the core spins 
are aligned will also have a fully polarized conduction 
band; therefore the crucial question becomes: what is 
the value of J? 

From the earliest days 1,3 through the revival of inter- 
est in the mid 1990s 9 and including recent literature 10-12 , 
many workers have assumed that J is much larger than 
the electronic bandwidth. Such a large value of J has 
been argued to lead to interesting observable effects in- 



cluding a temperature dependent optical spectral weight 
in the physically relevant lo < AeV regime 9 and to provide 
a direct connection between optical conductivity and the 
values of the magnetic transition temperatures 13 . A tem- 
perature dependent optical spectral weight has indeed 
been observed 14 . Even more importantly, the experimen- 
tal phenomenology 14 suggests that there is a strong con- 
nection between electronic energetics and magnetic order 
and that in particular changes in electron kinetic energy 
occurring as the system is driven through ferromagnetic- 
non-ferromagnetic transitions lead to the observed 15 mul- 
tiphase behavior, suggesting again a very large value of 
J. On the other hand, Measurements on gas-phase Mn 16 
suggest a relatively small J (comparable to the elec- 
tronic bandwidth), and measurements on related tran- 
sition metal oxides 17 ' 18 suggest that gas-phase values for 
J are only weakly renormalized by solid-state effects. 
On the other hand, many workers (including most of 
the references listed above) interpret experiments in the 
'double-exchange' picture which requires a J much larger 
than the bandwidth. Determining the magnitude of J 
and the effect of magnetic order on energetics is there- 
fore an important issue, which remains controversial. 

Tokura and co-workers noted 19 that the optical con- 
ductivity may contain direct information about the value 
of J: in a paramagnetic state optical processes were pos- 
sible in which an electron moves from one site, where it 
was locally parallel to the core spin, to an adjacent site in 
which (because the core spin is differently oriented) it is 
locally antiparallel. Such transitions should occur at an 
energy set by the Hund's coupling J, and should become 
visible as the temperature is raised. Ref. 13 analysed this 
issue in detail, presenting a comparison between theory 
and experiment which was argued to exclude the possibil- 
ity that the 'antiparallel' transitions had been observed. 

In this paper we reexamine the previous work 9 ' 13 , us- 
ing a more realistic tight- binding parametrization of the 
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relevant portions of the band structure in combination 
with a dynamical mean field treatment of the interaction 
between electrons and core spins and employing more 
detailed and extensive numerics. While various aspects 
of this physics have been extensively studied, especially 
for model systems, we believe that this paper is the 
first to combine a realistic band parametrization with 
a complete treatment of dynamical and thermodynamic 
quantities such as the conductivity (for both low-T and 
T > T c ), spin- wave stiffness and transition temperature. 
(Of course, many important studies of static properties 
using realistic bands have appeared: for a recent example 
see 20 . We also note that in early important work Taka- 
hashi and Shiba computed T — optical conductivities 
using a similar parametrization of band theory 21 .) We 
compute magnetic transition temperatures, optical con- 
ductivities and optical spectral weights and the T = 
spin-wave stiffness for a range of model parameters and 
by comparison to experiment estimate the actual phys- 
ical parameters for the material. We find that contrary 
to statements made in some previous work 13 a moderate 
value of J provides a reasonable account of many aspects 
of the optical spectrum and a plausible account of the 
magnetic transition temperature and T = spin wave 
stiffness. The differences with previous work apparently 
arise because the calculations of Ref 9,13 were performed 
on perhaps oversimplified model systems. Also some of 
the results 13 contain numerical errors. 



II. MODEL 



A. Hamiltonian 



We consider electrons moving among sites of a sim- 
ple cubic lattice according to a tight binding model with 
hopping matrix elements obtained from band theory and 
coupled to a local classical 'core spin'. Thus we write: 



H = H band + H int 



(1) 



For the interaction we take a local coupling of core 
spins to conduction electrons: 
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We obtain H^nd from LSD A band calculations 5,22 , 
which are consistent with straightforward quantum 
chemical (Goodenough-Kanamori) considerations. Both 
approaches suggest that the itinerant electrons lie in two 
symmetry related bands which may be thought of as 
arising from electrons hopping among e g symmetry local 
orbitals and are reasonably well described by a nearest 
neighbor tight binding model 13 ' 23 . To write this hopping 
it is convenient to adopt a Pauli matrix notation in which 
the up state is the \x 2 — y 2 > orbital and the down state 



is the \3z 2 — r 2 > orbital. The basic hopping is orbital 
dependent so that 



H 



band 



T E (^) a ' 



C- 



+*,a,« C *« + H - C -) ( 3 ) 



iSaba 



where S is a vector connecting site i to a nearest neigh- 
bor and the hopping is determined by a 2x2 matrix given 
after Fourier transformation by 



£ = £o(p) + ~e(p) ■ ~r (4) 
with ~t the usual Pauli matrices, 6 the lattice constant, 

£o(p) = -t(cos(p x b) + cos(p y b) + cos(p z b)) (5) 
and ~e(p) = (e x (p), 0, e z (p)) with 



-^(cos(p x b) - cos(p y b)) 
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The energy eigenvalues are 



(6) 



t(cos( Pz b) - -(cos(p x 6) + cos(p y b)) (7) 



E± = e ± \je 2 x + e 2 z 



(8) 



Note that along the zone diagonals ((1,1,1) etc) e z . x = 
so the two bands are degenerate and along the line 
to any cube face ((1,0,0) and equivalent) one of the two 
bands is dispersionless. These two features occur, to a 
high degree of accuracy, in the calculated band struc- 
ture, lending support to the nearest neighbor tight bind- 
ing modelling. 

The best fit to the calculated band theory leads to 
t = 0.67eV. This, for example, implies that at band 
filling n = 0.7 the density of states per Mn-ion is 0.57 /eV 
which reproduces almost exactly the local spin density 
approximation value 0.58 /eV per Mn-ion quoted on p. 
1154 of Ref 5 . (Note also that 5 contains a misprint 24 in 
the value of the Drude plasma frequency Q p = 1.9eV. 
The correct band theory value is very close to the £l p = 
y / A-Ke 2 D/b = 3. leV^ which follows using t — 0.67eV in 
Eq. 12 below). 

The energy bands extend from — 3< to +3t. Useful 
quantities to characterize the state of the system include 
the particle density 
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the density of states 
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FIG. 1: Upper panel: T = density of states for a ferro- 
magnetic ground state, plotted against chemical potential for 
several different Hund's couplings J. Lower panel: Density 
of states in paramagnetic phase plotted against energy for 
several J 



(whose T = value is shown for different J values in 
the upper panel of Fig. 1), the 'kinetic energy' per site 
(N is the number of sites) given by 



K _ <H band > 
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K= ^nH ta "W < 4a^i+Saa + H.C. > (11) 



and the 'Drude weight' D which for a spin polarized sys- 
tem and current flowing in the z direction is 
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The expectation value of the term in J proportional to 
A gives the total oscillator strength, S(po) in the con- 
duction band contribution to the optical conductivity 
( see i4,25, 26 for detai i s ). We have 



5(oo) = -K 



The conductivity is 



S{oo) - Xjjjfy 
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with Xjj the usual Kubo formula current-current corre- 
lation function evaluated using J evaluated at A = 0. 



C. Spin Wave stiffness 

To compute the T — spin wave stiffness we follow the 
standard procedure outlined in 27 ' 28 We compute the en- 
ergy cost of a small amplitude, long wavelength rotation 
of the order parameter away from the fully polarized fer- 
romagnetic state, which we take to be aligned with the z 
axis. The calculation is most easily carried out by locally 
rotating the spin quantization axis to align with the local 
spin direction, so that the T — Hamiltonian becomes 
(note repeated orbital and spin indices a, b, a, a, a' are 
summed over) 



H = -\ E t ah {5)ct +s ^c iM ,Rt +Saa R ia ^ + H.c. (16) 

i,6 

where R are the usual 5=1/2 rotation matrices. We 
find two terms. The first one arises from the term in 



Figure 2 shows the kinetic energy (divided by 3 to fa- 
cilitate comparison to optics) and Drude weight, plotted 
against particle density for a fully spin polarized gound 
state. 



B. Conductivity 

To obtain the conductivity we require a coupling be- 
tween the electric field and the electronic states. We rep- 
resent the electric field by a vector potential A and adopt 

the Peierls phase approximation, fj_j — ► ti-je 1 ? A ' Rij . 
This approximation has been argued to be accurate in 
other transition metal oxide contexts 23,25 and was also 
used by Takahashi and Shiba 21 . The current density 
operator in the z direction is 
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FIG. 2: Solid line: Kinetic energy (from Eq 11 but divided by 
3) in units of hopping parameter t plotted vs. particle density 
n. Dashed line: Drude weight D z /t (from Eq 12) 



R + R proportional to the square of the deviation, m q of 
the magnetization from its ordered state value and is (K 
is the kinetic energy defined above in Eq 11). 
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while the second one arises from inserting the expres- 
sion for the term in R + R which is linear in m q into the fa- 
miliar second order perturbation theory expression, and 
is 



(qa) m q ■ m,- q I S w(J) 



(18) 



with 



Isw(J) =£> 2sin2 (P*H*+(P) + *-(P)) (19) 



where (note we have suppressed the momentum labels to 
avoid clutter in the equations) 
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and tan(0) = e x j {\j E 2 . + e 2 + e^j . 

The spin wave stiffness Dsw as conventionally defined 
is given by twice the total coefficient of (q) 2 m q ■ m_ g 
divided by the ordered moment, i.e. 



Dsw = 



K 
12 



-2I SW (J) 



(22) 



This expression applies to the tight binding band struc- 
ture and to all J such that the ground state is fully polar- 
ized. It is consistent with those obtained by Kubo and 
Ohata 29 and Furukawa 30 who considered simpler mod- 
els in the infinite J and (Furukawa) first 1/J correction. 
We have evaluated D (Eq 22) for the chemical potential 
/i = — 1 corresponding to the widely studied 0.7 doping; 
results, made dimensionless by dividing by tb 2 are shown 
in Fig. 3. 



III. METHOD OF EVALUATION 

A. Overview 

To evaluate the non-ground-state properties of H 
we use the dynamical mean field method 13, 31-33 This 
method is extensively described and justified elsewhere, 
and is relevant here because the principal interactions 
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FIG. 3: Dimensionless spin wave stiffness Dsw /(tb 2 ) as func- 
tion of Hund's coupling J divided by bandwidth t, for T = 
fully polarized state of two-orbital lattice model at carrier 
concentration n = 0.7, calculated using Eq. 22 



are local. In brief the central approximation concerns 
the electron Green function. For the band structure of 
present relevance the electron Green function is written 
in general as 



G(z,p) = {z-S(p,z)-e(p) + ^y 



(23) 



with S(p, z) the self energy. In the dynamical mean 
field method one takes S to be p— independent, i.e. to 
depend only on z. The important quantity is then the 
momentum integrated Green function 
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mom—int 



(z) = b 3 J 



d 3 p 



G(z,p) 



(24) 



The local Green function, being a function of frequency 
only, may be derived from a local quantum field theory, 
which is specified by a partition function Zi oc given in 
terms of a mean field function a(r). This is in general 
a matrix in orbital and spin indices (which we do not 
explicitly write here) 



Z loc = J Vc+cExp[J drdr' c + (t)&{t - t')c(t') 

+ J drH mt ] (25) 

^From this action one may extract a local Green func- 
tion Gioc and self energy S via 



G 1oc (t) = 



SlnZioc 
Sa(r) 



(a-FJ)- 1 



(26) 



The mean field function a is fixed via the constraint 
that the local Green function calculated from Zi oc is 
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identical to the momentum integrated Green function ob- 
tained from Eq 24, using the local self energy defined in 
the second equality in Eq 26, i.e. 



(27) 



In the present problem the two orbitals are degenerate 
so it is not necessary to consider orbital indices in G\ oc . 
For T=0K a ferromagnetic core spin configuration yields 
a spin dependent but frequency independent S CT = <j z J. 
For T > T c Eq 27 is an integral equation for S, which is 
solved by numerical iteration. 



B. Conductivity 



In the dynamical mean field approximation there are 
no vertex corrections for the current current correlation 
function 33,34 , so it is given by 



d 3 P 



Tr[3 z G(p,iQ + iuj n )3 z G(p,iuj n )} (28) 

with G given by Eq. 23 using the momentum indepen- 
dent self energy from the final solution of Eq 26. 

The Matsubara sum and analytical continuation may 
easily be performed, and the dissipative part of the con- 
ductivity is given by a(Q) = Im[xjj(Sl — iS)]/Q. It is 
convenient to transform within the trace to a orbital ba- 
sis which diagonalizes the Green functions. In the model 
considered here the self energy is the same for both or- 
bitals, and so we may split the result into an intra- and 
interband part, writing 
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FIG. 4: Optical conductivity for the ferromagnetic ground 
state at n=0.7, computed as described in the text using 
t — 0.67eV and a single-particle scattering rate 1/r = O.leV. 
Note the change in y-axis scale at leV. 



in the equations: for example, the band-n (n = 1,2) 
spectral function A ( - n \p,u>) = Jto[E(u;)]((u; — E n (p) — 
i?e[S(w)]) 2 + /m[S(w)] 2 )- 1 . 

At T > T c the conductivity is calculated using the 
self energy obtained from the solution of the mean field 
equations. For T — > the spin-induced self energy van- 
ishes (except for the spin-dependent energy shift) and 
we have added a modest impurity broadening (scattering 
rate S"(w) = r _1 = O.lt) so the Drude peak is visible 
on the same scale as the interband term. The calculated 
T = conductivity is shown in Fig. 4. 
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The quantities A, e and $ depend on momentum p which 
we have not written explicitly to avoid confusing clutter 



C. Transition Temperature 

To determine the transition temperature it is conve- 
nient to write the mean field function in the presence 
of a uniform magnetization m, which is in principle a 
matrix in spin space, as 



a = a l + a\o ■ m 



(33) 



and to take the magnetization direction to be parallel 
to z so that after integrating out the fcrmions we ob- 
tain for the impurity model action (the two is for orbital 
degeneracy) 



6 



Simp = 2j2 Tr M fl o + a\<j z + J~S ■ ~a] (34) 

n 

We can evaluate the trace over spins, getting (9 is the 
angle between the core spin direction and the magneti- 
zation direction, which is taken to be the z direction) 

S imp (9) = 2 M a o - a \ - j2 S 2 - 2a i JS C0S W] ( 35 ) 

n 

Assuming (as detailed studies of simpler model systems 
have confirmed 9 ) that the transition is second order we 
may study it by linearizing in a±, so that the magnetiza- 
tion m is given by 



fd cos 9 {cos(6)Exp[S '(0)]) 
m ~ < J d cos 9 Exp[S (9)} > 

+ 0(4) 
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FIG. 5: Calculated magnetic transition temperature plotted 
against band filling for several different J values, along with 
J = oo Bethe-lattice result (for comparison to previous work). 



The impurity (local) Green function corresponding to 
Eq 35 is 



Gioc = ( (a + a\a z + J S ■ ) 

(a - cr z [»l(l + | a /j2 S 2 ) + Jm]) 
= («o - j" 2 S 2 ) 



+ 0{a\) 
(37) 



where the average has been taken over the action defined 
by Eq 35. From Eq 37 we may compute the self energy 
S = -Gf l + a + ai<7 z finding S = S + Sicr z : 



J 2 S 2 

a 
1J 2 S 2 



So 
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Finally, a is fixed from 



Gloc = \Tr ( — ^ [w - S - SiC7 z - e p ] 1 
2 7 (27r) d 

1 f d d p , ,_i 

[w - S - e p ] _1 Sio-2 [w - S - e p ] _1 
= /i(w-So) + S lf 7 z / 2 (w-So) (40) 

-th /„(,) = J ^(^)". 

Eqs. 36-40 give an equation for T c , which contains 
only ao and the paramagnetic solution of Eq. 27 which 
we denote here as Go- 



T C = -AJ 2 T C J2 



Gl-h 



(3a 2 - J 2 )G 2 + J 2 h 



(41) 



This formula shows explicitly that in weak coupling 
T r ocJ 2 . 



IV. RESULTS 

As noted above, the band calculation is well fit by a t 
of about 0.67eV. The physically relevant densities corre- 
spond to n < 1 and the most experimentally relevant den- 
sity range is n — 0.6 — 0.8 where the ground state is in fact 
a ferromagnetic metal. We have calculated the magnetic 
transition temperature as a function of band filling for 
different values of the Hund's coupling J. Representative 
results are shown in Fig. 5. We see that the calculated 
transition temperatures become noticeably higher than 
the experimental range (~ 400K) once J becomes greater 
than a number of the order of unity. Ref 13 argued that 
the calculated transition temperatures for J >• t were 
consistent with experimental data. We now believe that 
this conclusion was based on normalizing to an incorrect 
value for the kinetic energy and should be disregarded. 
The present calculation reveals, in agreement with results 
obtained by previous workers, that even modest values 
of J suffice to push transition temperatures into a range 
higher than the highest observed T c . Corrections to the 
mean field approximation arc believed 13 in this problem 
to reduce T c by of the order of 30% and do not change 
this conclusion. 

The transition temperature data therefore suggest that 
a relatively modest J, corresponding to a local spin flip 
energy 2J ~ 2 — 3eV and less than the full bandwidth 
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FIG. 6: Optical conductivity in the paramagnetic phase at 
(T > T c ) for different couplings J and density n=0.7, com- 
puted using t = 0.67eV. 
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FIG. 7: Main panel: difference a(u,T > T c )-o(u>,T = 0) for 
n = 0.7 and J = 0.67eV (solid line) 1.33eV (dotted line) J = 
2eV (long dashed line), J = 2.67eV (second solid line) and 
J = 3.33eV (dot-dashed line), inset: difference conductivity 
as function of density for J = 1.33eV (n = 0.5 (dotted), n = 
0.6 (solid) n = 0.7 (dashed) n = 0.8 (dash-dot) and n = 0.9 
(second solid line). All curves computed using t = 0.67eV. 



~ 4eV, might be appropriate for the CMR manganite 
materials. Such a value of J would be consistent with 
the value J w 2.7eV appropriate for gas-phase Mn (note 
that spectroscopic studies such as 17 ' 18 indicate that the 
Hund's coupling is much less renormalized by solid-state 
effects than is the on-site Coulomb interaction). Our 
computed values for the T = spin wave stiffness also 
support a lower value for J. Values in physical units may 
be obtained by multiplying the data shown in Fig 3 by 
tb 2 9.67 (using the numerical values t — 0.67eV and 
a = 3.81). A J — 2.5t yields a D = 300meVA 2 , already 
much higher than the D 160 — 190eVA 2 observed 35 in 
optimally metallic manganites. 

Further evidence that a smaller value of J might be 
appropriate comes from our calculated optical conduc- 
tivity, shown in Fig. 6 for T > T c and various values of 
J. Comparison to Fig 4 shows that if J is larger than 
about leV, then as T is raised a new feature becomes 
visible at an energy of order 2 J. This feature arises from 
motion of carriers from one site to another site with an 
antialigncd spin, and its existence was noted by Tokura 
and co-workers 19 . Such a feature is visible in experimen- 
tal data at an energy of about 3eV 14 . Comparison of 
Fig. 7 to the data of 14 shows that a J = 2.5t w 1.5eV 



would approximately reproduce the observed peak posi- 
tion and magnitude. For smaller values of J the peak 
becomes difficult to distinguish from the interband tran- 
sition; for larger values the peak moves to higher ener- 
gies than is observed and the oscillator strength decreases 
below what is observed. The Table shows the oscillator 
strength in the " peak at 2J", obtained by integrating 
<7t c — Co from the highest zero-crossing up to w = oo 
and converting to kinetic energy via Eq 14 (e.g. the in- 
tegration for J = 3.33eV and n = 0.7 includes the area 
between 5.8eV and 9.6eV). ft is also convenient to define 
Kt c , the total spectral weight of a at T c , by evaluat- 
ing Eqs 11, 14 at T = T c . The ratio K anti /K Tc can be 
compared with the experimental results of Rcf 14 . 



Ref argued, on the basis of calculations based on an 
oversimplified Bethe-lattice model, that the observed os- 
cillator strength in the putative 'peak at 2J' was too 
large to be consistent with the calculation. The present 
calculations, which are based on the more realistic tight 
binding model, suggest that this is not the case. 

Previous work 9,13 has drawn attention to a relation- 
ship between optical spectal weight (i.e. kinetic energy) 
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FIG. 8: Comparison of optical spectral weight (expressed as 
a kinetic energy) at Tc (dashed and dotted lines) and K 
(solid line) and various J; the minimal J value, which gives a 
halfmetallic ground state is dependent on n (e.g. for n=0.7 is 
J mjn =0.64eV) 



changes and the value of the magnetic transition tem- 
perature. Our new results indicate that this result is of 
less generality than suspected. Fig. 8 compares, for dif- 
ferent values of J, the optical spectral weight at T = T c 
(dashed and dotted lines) and T = (solid line). The 
J = oo limit reveals the expected approximately 30% de- 
crease in spectral weight between T = and T = T c , but 
the other traces reveal that the situation is more compli- 
cated in general. Indeed, at T — and for J sufficiently 
large, only one spin band is occupied, and the band is 
therefore filled to a higher level, thereby losing kinetic 
energy. On the other hand, at T > T c both spin direc- 
tions are equally probable, so electrons can redistribute 
themselves to the lower parts of the band, thereby gain- 
ing kinetic energy. However, in the paramagnetic phase 
the non-vanishing J blocks some hopping processes. At 
intermediate J these two processes compete, and as seen 
in Fig. 8, which one is dominant depends on both the 
value of J and the band-filling. The blocking process is 
always dominant only at very large J. 

The non-systematic change of kinetic energy upon 
magnetic ordering observed for reasonable values of J 
suggests that the relation, proposed in previous work 13 , 



TABLE I: the kinetic energy (oc spectral weight), which corre- 
sponds to the "spin flip" transition in the paramagnetic phase 
presented in absolute units and relative to the total T = T c 
spectral weight 



J[eV] 


—Kanti[meV] 


Kanti/ Ktc 


0.67 


13.3 


0.038 


1.33 


55.6 


0.17 


2.00 


56.7 


0.19 


2.67 


48.9 


0.17 


3.33 


42.2 


0.15 




-0.05 0.05 0.1 0.15 0.2 0.25 

K(T=0) - K(Tc) [ eV ] 



FIG. 9: Upper panel: Open dots: calculated transition tem- 
perature plotted against change in kinetic energy between 
T — T c and T = for different J as shown with density 
n as implicit parameter. Dashed lines are guides to the eye. 
Along each line the doping increases by 0.1 steps along the 
lines from 0.1 up to 1 or 0.9. Lower panel: Filled dots: calcu- 
lated transition temperature plotted against change in kinetic 
energy between T — T c and T — for different fixed densi- 
ties, as shown, with Hund's coupling J as implicit parameter. 
Dashed lines are results of a linear regression for the T c vs 
AK(J) curves. Solid line: results published in 13 for a Bethe 
lattice; the difference seems to arise because of an omitted 
factor of two in 13 . 



between change in kinetic energy and value of transition 
temperature is not as general as expected. The lower 
panel of Fig. 9 shows the value of the transition temper- 
ature, plotted against change in kinetic energy for several 
different carrier densities, with J as an implicit parame- 
ter. Also shown, as the solid line, is the result obtained 
in Ref 13 (the different slope occurs because the results 
of 13 are for an orbitally non-degenerate model). A lin- 
ear dependence seems to be reasonable and is consistent 
with the ideas in Ref. 13 - but we observe, that the slope 
is doping dependent. The upper panel in Fig. 9 shows 
the same plot for several different J, with carrier density 
as an implicit parameter. A much more complicated ap- 
pearance is obvious. We see that the relation between 
T c and kinetic energy is not at all universal, but depends 
strongly on electron density, essentially for the reasons 
given in the paragraph above. 
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V. CONCLUSIONS 

We have presented a detailed analysis of the opti- 
cal conductivity and magnetic transition temperature 
of electrons using a quasirealistic band structure and 
Hund's coupling to classical 'core spins'. Our main find- 
ing is that, in contrast to a widely-made assumption in 
the manganite literature but in agreement with quantum 
chemical estimates and gas-phase measurements, a mod- 
erate value of the Hund's coupling J w 1.3— 1.5eV = 2.5t 
(large enough to fully polarize the conduction band at 
T = but less than the LSDA band width) suffices to 
account both for the order of magnitude of the magnetic 
transition temperature and for a feature observed in the 
optical conductivity in the paramagnetic phase. We have 
also presented predictions for the evolution of this feature 
as the carrier density is varied. Errors in previous work 13 
(arising from use of an oversimplified band structure and 
from numerical errors) have been corrected. We have also 
presented results for the spin-wave stiffness expected for 
a fully polarized ground state and a tight-binding band 
structure at arbitrary J. 

We observe that even with a modest J the spin wave 
stiffness and transition temperature are rather overesti- 
mated. It is possible to fix up these discrepancies by 
adding a phcnomenological antifcrromagnetic core-spin- 
core-spin coupling, as has been proposed by many au- 
thors. However, in our view a more serious problem re- 
mains: if these estimates for J are accepted, then our ba- 
sic understanding of the physics of the manganites must 
be reconsidered, because for these value of J the change 
in kinetic energy between ferromagnetic and paramag- 
netic state is very small. A series of experiments has 
made it clear that the 'colossal' magnetoresistance is as- 
sociated with a change in electronic state from a more 
or less fermi-liquid like state to a state with strong lo- 
cal lattice distortions, leading to either short ranged or 
long ranged charge and orbital ordering. Previous work 9 , 
based on a large- J limit, had predicted this behavior as 
a consequence of the decrease in kinetic energy as spins 



were disordered, allowing the localizing tendency of the 
electron-phonon interaction to overcome the banding ten- 
dency of the electrons and produce a new phase. An 
approximately 40 per cent decrease in conduction band 
oscillator strength, which is naturally associated to a de- 
crease in electronic kinetic energy, is indeed observed, 
and was argued to be consistent with this physics. How- 
ever, if J is in the range proposed here, then the dis- 
ordering of the spins does not significantly change the 
electronic kinetic energy and another explanation for the 
high temperature second phase and thus for the 'colossal' 
magnetoresistance must be sought. 

One possibility is that a 'Hubbard U' interaction is 
important, and reduces the electronic kinetic energy and 
band width to the point where even the modest J we find 
here can have a strong effect. However, the total con- 
duction band spectral weight observed at low T is close 
to the full band theory value, and indeed the while our 
simple calculation does not reproduce the detailed line- 
shape, the calculated magnitude of the conductivity in 
the l-3eV is in reasonable agreement with data, sugg est- 
ing that Hubbard-U effects are perhaps not so important. 
It is possible, however, that the good agreement is acci- 
dental, and that contributions from other orbitals (for 
example, 'wrong spin' t2 g electrons) are important in the 
real materials, although they are absent in our calcula- 
tions. We also note that the susceptibility of the finite J 
paramagnetic state to local (or long ranged) lattice dis- 
tortions has not been computed. It might be that even 
though the kinetic energies of the T = and T > T c 
phases are similar (at moderate J), the T > T c phase is 
more susceptible to local lattice distortions. This possi- 
bility urgently deserves further investigation. 
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